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Abstract 

We investigate the effect of orbital overlap on optical matrix elements in empirical tight-binding 
models. Empirical tight-binding models assume an orthogonal basis of (atomiclike) states and a 
diagonal coordinate operator which neglects the intra-atomic part. It is shown that, starting with 
an atomic basis which is not orthogonal, the orthogonalization process induces intra-atomic matrix 
elements of the coordinate operator and extends the range of the effective Hamiltonian. We analyze 
simple tight-binding models and show that non-orthogonality plays an important role in optical 
matrix elements. In addition, the procedure gives formal justification to the nearest-neighbor spin- 
orbit interaction introduced by Boykin [Phys. Rev B 57, 1620 (1998)] in order to describe the 
Dresselahaus term which is neglected in empirical tight-binding models. 

PACS numbers: 73.21.Cd, 73.21.Fg,78.64.De 
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I. INTRODUCTION 



The tight-binding (TB) approach to electronic structure is one of the most used methods 
in sohd state systems^. The empirical tight-binding (ETB) method, which dates back to 
the work of Slater and Koster-, assumes mostly two-center approximation and the matrix 
elements of the Hamiltonian between orthogonal and atom-centered orbitals^ are treated 
as parameters fitted to experiment or first-principles calculations. ETB is widely employed 
to the description of electronic structure of complex systems^ like interfaces and defects in 
crystals, amorphous materials, nanoclusters, and quantum dots because it is computationally 
efficient (up to 3 order of magnitude faster than the ab intio density functional methods) 
and provides physically transparent results. Many calculations consider just the nearest- 
neighbor Hamiltonian with fewer parameters but with additional orbitals introduced^. To 
consider a higher accuracy, the range of the Hamiltonian is extended to few nearest-neighbor 
shells (up to the first three shells) and, therefore more fitting parameters^. However the use 
of a non-orthogonal formalism might scale back the range of the Hamiltonian having the 
additional fitting from the overlap matrix. In some instances^, when strain is present, non- 
orthogonality is invoked implicitly to accommodate the changes of the on-site energies due 
to local displacements in addition to the well known scaling of the transfer integrals^. A 
non-orthogonal formalism has also a less obvious advantage. Because they have a longer 
range than the atomic orbitals, the orthogonalized orbitals samples the local environment, 
making them better suited for transferability to complex systems^. 

Calculation of optical spectra in the ETB formalism requires the knowledge of additional 
parameters: the momentum or velocity matrix elements between initial and final states. In 
the early work, momentum matrix elements were considered as extra parameters fitted to 
the experimental or first-principles calculated dielectric function. However, ETB has been 
extended to include the interaction with electromagnetic fields^'^ by making the substitution 
p = C^o/f^) VkH, such that dielectric function and other optical properties can be calculated 
without additional parameters. The scheme is based on the Peierls substitution of Hamilto- 
nian matrix elements^i allowing us to calculate directly the momentum or velocity matrix 
elements. In Ref. Q andQ it is shown that the substitution p = i^^/fi) VkH leads to 
the neglect of the intra-atomic momentum matrix elements or, equivalently, the coordinate 
operator is diagonal in the subsequent basis as we will indicate below. However, the Peierls- 
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tight-binding (i.e. zero intra-atomic position parameters) has been successfully used in Ref. 
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Pedersen et al^ introduced an additional momentum matrix element to accommodate 
the intra-atomic transitions. In contrast, Boykin and VogU^ showed that adding intra-atomic 
terms suppresses the gauge invariance. To circumvent this problem B.A. Foreman^^ used 
group theory arguments to construct the basis in which intra-atomic matrix elements are 
present and the lattice gauge theory to define the interaction of electromagnetic fields with 
electrons in crystals. 

The effect of orbital overlapping on electronic structure has been studied for simple 
system a^^i^^ . In this paper we investigate the optical matrix elements in the presence of 
non-orthogonal (overlapping) orbitals. As far as we know, no study has been done in 
this direction. We show that intra-atomic contributions of the coordinate operator are 
induced simply by the orthogonalization process. The orthogonalization process induces 
terms equivalent with more distant interactions, such that it gives formal justification for 
the nearest-neighbor spin-orbit interaction introduced in Ref. IJJI for TB model with spin- 
orbit interaction^ii. The analysis of simple systems shows that the non-orthogonal orbitals 
play an important role on optical matrix elements. We re-analyze the example of Pedersen 
et al}^ to show that the non-orthogonal orbitals improve the optical matrix elements. In 
the case of graphene, the overlap and TB parameterization are crucial in explaining the ex- 
perimental data. Moreover, similar arguments can be employed in the ab-initio TB-LMTO 
(tight binding linear muffin-tin orbitals) method^i, leading to faster calculations of optical 
matrix elements in a parameter free theory. 

II. TIGHT-BINDING CALCULATIONS AND NON-ORTHOGONALITY 

To fix ideas we consider a localized basis | aR) , where a is the orbital type and R is the 
center of the orbital (Lowdin orbitals)^. The crystal Hamiltonian H is diagonalized within 
the Bloch sums of the localized basis 

|afc) =y^J]e^^^|«/?) (1) 

R 



3 



as follows 



\nk) = ^Cnaik)\ak), (2) 



with 



H {k)\nk) = Enk\nk) . (3) 
The kinematic momentum operator involved in optical transitions is defined as 

P='^[r,H]- (4) 
In the crystal momentum representations?,^ the kinematic momentum operator is 

ITl 

P=jVkH{k), (5) 

where H (k) is the Hamiltonian in the crystal momentum representation. Eq. (0) holds in 
a complete basis as well as in an incomplete basis. However, in an incomplete basis the 
momentum operator p and coordinate operator r do not satisfy the canonical commutation 
relations leading to different formula for effective masses and Peierls-coupling formula in- 
volving the vector potential. These issues are detailed in Refs. 
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23 and^. The coordinate 



operator r is considered to have the following matrix elements in the localized basis | aR) 



{a'E! \r\ aR) = {R5aa' + d^a') Srr', (6) 

since the overlapping of the orbitals belonging to different atoms is supposed to be small. 
Here d^a' is the intra-atomic matrix element. In the usual tight-binding theory the coordi- 
nate operator is diagonal^''. Therefore the intra-atomic parts are neglecte d^^i^^ leading to no 
need of other fitting parameters beyond those of the Hamiltonian and to gauge invariance. 

Pedersen et alM pointed out that there are cases in which the neglect of the intra-atomic 
part may conduct to the underestimation of the momentum operator arguing that by using 
Eqs. dH- (jlD, 

{nk\p \mk) = '-^ ^ c w (k) Cma (k) Vfc {a'k \H\ ak) + 

(7) 

~ {^nk S'mfc} C *na' {k) Cma (^) da'a 

a,a' 
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The neglect of the second term in Eq. ((Tj) reproduces Eq. (0). Therefore in the tight- 
binding basis, which is finite, by using Eq. (jS)) one is neglecting the second term in Eq. ((7j) 
or the intra-atomic partf^^. This shortcoming happens because the momentum and position 
operators do not satisfy the canonical commutation relations in a finite basia^i^. One way 
to add intra-atomic terms is the construction of B. A. Foreman^^. However, intra-atomic 
terms can be induced if one considers a non-orthogonal basis. To show this, let us have an 
atomic basis with non-zero overlapping 

(XaR I Xa'R') = 1 + Saa'RR' = 1 + S. (8) 

The orthogonal basis corresponding to Eq. (jH)) is (the Lowdin procedure) 

Ix') =(1 + S)-'^|X) (9) 
In the new orthogonal basis an operator transforms according to 

A' = {1 + S)-^^A{1 + S)-^^ . (10) 
Formally, expanding Eq. in power series of S we rewrite Eq. fllOj) as 

A' = A-^{SA + AS) + ^ {ASS + SSA) + ^SAS .... (11) 
The inverse transform of Eq. (fTUI) has the following expansion 

A = A' + -{SA' + A'S) -- {A'SS + SSA') + ]sA'S ... (12) 
2 8 4 

Now suppose that in the Lowdin basis the intra-atomic matrix element daa' is zero, such 
that the Hamiltonian fulfils the gauge invariance conditions. In the original non-orthogonal 
(atomic) basis, however there are intra-atomic elements. These can be easily seen if one 
applies the inverse transform Eq. ()12|) (Fig. ^). Thus in the atomic basis, up to the second 
order in S, the intra-atomic matrix element is 



^RaRa' — ^ ^ {—rRSRa',R"a"SR"a",Ra — SRa',R"SR"a",RarR + 2Sjia'^R"a"rR"Sii»a",Ra) 



R" 



(13) 
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Eq. ()13|) shows us that in the atomic basis the intra-atomic matrix elements of coordi- 
nate operator are non-zero. Hence, considering the overlap, intra-atomic optical transitions 
can be incorporated. Although the intra-atomic corrections are second order in the over- 
lap S, the overall corrections to the optical matrix elements are first order in S. In the 
same time the range of the Hamiltonian has been increased by applying the transformation 
given by Eq. (fTT|) to the Hamiltonian matrix (Fig. ^). This result suggests that although a 
nearest-neighbor Hamiltonian might give a good reproducibility of the electronic structure, 
it completely misses the intra-atomic terms of optical matrix elements. The relationship 
between the overlap and longer ranged Hamiltonians is able to explain the nearest-neighbor 

n 

spin-orbit interaction introduced in Ref. 1191 in order to reproduce the Dresselhaus terms 
in zinc blend structures. Thus, the spin orbit contribution to the optical matrix elements, 
VkHso {k), is non-zero. In the same time the initial prescription given by Chadi^S is pre- 
served. The overlap and long-range Hamiltonians are also closely interrelated in quantum 
wire transport. Using nearest-neighbor Hamiltonians, the overlap is crucial in explaining 
anti- resonances in quantum wires^^. However the same anti-resonances are reproduced with 
a Hamiltonian in which the effect of overlapping has been transferred to the second nearest- 
neighbor hoping elementa^^. 

In a recent paper— it is shown that a piece- wise constant coordinate operator (and there- 
fore diagonal) in ab-initio TB-LMTO (tight binding linear muffin-tin orbitals) methods is 
analogous to the coordinate operator in semi-empirical methods. The most localized rep- 
resentation (TB representation), where the Hamiltonian is short-ranged, is not the best for 
calculations although it is advantageous for numerical treatments. On the contrary, the co- 
ordinate operator was considered diagonal in the (nearly) orthogonal representation (with a 
long-ranged Hamiltonian) and used in transport calculations. The results were in very good 
agreement with the experimental values and with the results with the exact evaluation of the 
coordinate operator. Thus non-orthogonality plays an important role not only in empirical 
models but also in first-principles methods. If one assumes that the coordinate operator is 
piecewise constant and that the assumption is good enough, the calculations of the opti- 
cal matrix elements can be obtained faster from electron band calculations. From Eq. ((Tj), 
the A;-derivative of the Hamiltonian is calculated by fast Fourier transformations. Thus it 
is more computationally efficient than the usual scheme presented in Ref. |27l However, 



the applicability of a piecewise constant coordinate to optical properties of various physical 
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FIG. 1: (a) Schematic representation of the intra-atomic matrix elements of the coordinate oper- 
ator induced by the orbital overlap; the generic model has two orbitals per site, s-like (-) orbitals 
and p-like (- -) orbitals . We illustrate the matrix elements r of the coordinate operator in the 
orthogonal basis by small dash arrows , the overlap matrix elements S by dash arrows, the matrix 
elements d of the coordinate operator in non-orthogonal basis by full arrows, and the Hamiltonian 
matrix elements H in the non-orthogonal basis by dotted arrows, (b) Schematic representation of 
the increase in the range of the Hamiltonian in the orthogonalized basis by using the first order ap- 
proximation in Eq. The new long-range matrix element of the Hamiltonian in the orthogonal 
basis connecting site with site 2 for a nearest-neighbor Hamiltonian is shown by full arrow. 



systems remains to be investigated. 



III. OPTICAL MATRIX ELEMENTS IN SIMPLE TIGHT-BINDING MODELS 
WITH OVERLAP 

In the following we analyze the one-dimensional monoatomic crystal with two orbitals 
per atom, the one- dimensional diatomic crystal with one orbital per atom, and the two- 
dimensional graphene. 

a. monoatomic chain with two orbitals per atom. Schematic representation of a 
monoatomic chain with two orbitals per atom is given in Fig. ^ In a Bloch basis con- 
structed from the overlapping orbitals, the nearest-neighbor tight-binding Hamiltonian for 
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FIG. 2: (a) Band structure of the exact Kroning-Penney model and its approximations with a 
tight-binding model with overlap (see the text). Notice that following the same fitting procedure, 
the approximate bands are indistinguishable, (b) Momentum matrix elements of the exact Kroning- 
Penney model and its approximations with a tight-binding model with overlap. Notice the large 
variation of the momentum matrix elements with respect to the overlap while the bands are almost 
identical. 

a monoatomic chain with two orbitals per site is a 2x2 matrix 



Hik) 



(14) 



Es + 2 Vss cos (kL) 2 iVsp sin (kL) 
-2zV5psin (kL) Ep + 2 Vpp cos (kL) 
where Es and Ep are the energies of s-hke and p-like orbitals, respectively, Vss and Vpp are 
the coupling of two nearest neighbor s-like and p-like orbitals, respectively, and Vsp is the 
coupling of a s-like orbital with the nearest neighbor p-\ike orbital. L is the length of the 
unit cell and k is the wave vector. The overlap matrix has a similar form 



Sik) 



(15) 



1 + 25*55 cos (kL) 2iSspsin (kL) 
—2iSsp sin (kL) 1 + 2Sppcos (kL) 

In the orthogonal basis constructed according to Eq.(ini), the Hamiltonian matrix H is given 
by Eq. (fTn|l . The electronic bands are given by solving the eigenvalue problem H (k) \ nk) = 
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Enk \nk) and the interband matrix element of the kinematic momentum operator is p{k) = 
^ H' (k) \2k) , with H' the derivative of H with respect to k. We apply the above 
model to approximate the lowest two bands of the one-dimensional Kronig-Penny model. 
The Kronig-Penney model is a set of quantum wells of width a separated by barriers of 
height Vo and width b. The case is investigated by Pedersen et al}^ to suggest the need 
for intra-atomic contributions to optical transitions. We consider their strong-coupling case 
with a =8 A, h =1 A, and Vq = b eV. The first state in the quantum well is an s-like state, 
while second state is a p-like state. Accordingly, in the tight-binding counterpart of the 
Kronig-Penney model, the overlap matrix elements Ssp and Spp have to be negative. We 
adopt the same procedure^ for fitting the energy bands of the Kronig-Penney model. The 
absolute values of the overlap matrix elements are chosen to be the same for Sssi Ssp, and 
Spp. The results are shown in Fig.|21for an overlap of 0, 0.03 and 0.05 in comparison with the 
exact results of the Kronig-Penney model. While the energy bands are indistinguishable for 
tight-binding counterparts and agree well with the exact values, the interband momentum 
matrix elements vary and move toward exact values of the Kronig-Penney model. Because 
the absorption spectra are determined by the square modulus of the momentum matrix 
elements the above result is quite remarkable in the following sense as we explain below. 
Although we considered the strong coupling case (thin barriers) , the coupling between s-like 
and p-like states is weak (the matrix element Vsp is an order of magnitude smaller than the 
other matrix elements) such that the electron bands have almost either s- or p-like character 
over the entire Brillouin zone. Therefore, Vsp determines the magnitude of the interband 
momentum matrix element. In the same time the validity of p{k) = ^ {lk\ H' (k) \2k) is 
appropriate for strong inter-atomic coupling, such that the nearest-neighbor tight-binding 
model with orthogonal orbitals is inappropriate to calculate optical properties for the above 
model. Finally we want to mention that in one- dimensional crystals with inversion symmetry 
the coordinate operator is diagonal in the basis generated by the Wannier functions^S,. Hence, 
the "closer" to the Wannier functions are the Lowdin orbitals, the better reproduced are 
the momentum matrix elements. 

b. One- dimensional diatomic crystal with one orbital per atom. The chain is repre- 
sented by s-like orbitals at positions nL and p-like orbitals at nL + L/2, where L is length 
of unit cell and n is integer. The interaction up to the second-nearest neighbor is illustrated 
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FIG. 3: Schematic representation of the diatomic Hnear chain with one orbital per atom and 
lattice constant L . The first type of atoms are depicted as small and full circles and the second 
type as stripped circles. The interactions between atoms are shown by arrows. 

in Fig. ini The corresponding Hamiltonian matrix is 



H 



(16) 



Es + 2Vss cos {kL) 2iVspsin {kL/2) 
-2iVspsin {kL/2) Ep + 2Vpp cos (kL) 
Similar form holds for the overlap matrix. This can be an approximate model for super- 
lattices of type II, such as InAs-GaSB. In the InAs-GaSB superlattice the central feature 
is that the top of the GaSb valence band lies higher in energy than the bottom of the 
InAs conduction band, such that the electron and hole wave functions are overlapping. The 
electron/hole wave function is modeled by s-like/p-like orbitals. Keeping only the nearest 
neighbor interaction and overlap, the effect of overlapping is to decrease the momentum 
matrix elements as it is shown in Fig. |3J This simple result might help in explaining the 
increase of the photoluminescence intensity with the reduction of the electron-hole wave 
function overlap^^, which is not explained by the empirical pseudo-potential calculations 
used to for this purpose^A. The empirical pseudo-potential method2S used is non- atomistic, 
i.e. in their approach the Hamiltonian of the InAs/GaSb superlattice is constructed from 
the potential form factors of the InAs and GaSb bulk constituents. The potentials of the two 
bulk constituents are matched continuously at the interfaces such that there are no In-Sb 



or Ga-As bonds at the interface as there must be. As pointed out in Ref. 
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FIG. 4: (a) Optical matrix elements of the tight-binding model with overlap for a diatomic linear 
chain with one orbital per atom, (b) Energy bands of the tight-binding model with overlap for 
a diatomic linear chain with one orbital per atom. Full line is for overlap, dash line for a 0.01 
overlap, dot line for a 0.05 overlap, and dash-dot line for a 0.1 overlap. We choose arbitrary units 
because the system is rather generic. 

description is desired to take into account charge redistribution, segregation, and interdiffu- 
sion at the interface between InAs and GaSb. In contrast to Ref. I^O, Magri and Zungei^ 
solve the single-particle Schrodinger equation for each atom in the structure making their 
method atomistic. In this sense, TB models preserve the atomistic description of interfaces. 

c. Graphene. Recently, graphene as a two-dimensional sheet of graphite has been 
widely studied in the context of carbon nanotubes^. Graphite consists of a stack of graphene 
sheets , piled up and weakly interacting one with each other. Graphene has a hexagonal 
structure with two atoms in the unit cell (Fig. and very strong sp"^ bonds, causing a 
threefold coordinated planar structure. The remaining p^ orbitals are perpendicular to the 
plane, forming tt (bonding) and tt * (antibonding) states. The overlap of it electrons with the 
intra-plane sp"^ orbitals is small and vr and tt * electronic states dominate the physical prop- 
erties at low energy, around Fermi level. From Fig. we easily deduce the nearest-neighbor 
tight-binding Hamiltonian and overlap matrix for tt and tt * states as 
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FIG. 5: (a) The unit cell of two-dimensional graphite is shown as the dotted rhombus, (b) The 
Brillouin zone of two-dimensional graphite is shown as the shaded hexagon, aj, and bj, (i =1;2) are 
basis vectors and reciprocal lattice vectors, respectively. Energy dispersion relations and optical 
matrix elements are calculated along the perimeter of the dotted triangle connecting the high 
symmetry points, T, K, and M. 



^ To/- (k) Ep 

and 



Sik)^\ \/°^(^)|, (18) 

sof* (k) 1 

with / (A;) = e''"""/"^ + 2e"'''''"/^'^cos (^), 7o is the nearest neighbor transfer integral, Ep 
is the energy of tc orbitals, sq is the nearest neighbor overlap integral, a(=0.246 nm) is the 
lattice constant of graphite, and k is the two dimensional wave vector. Experimental data or 
first principles calculations put 70 between 2.5 and 3 eV, Ep = eV , and sq is found to be 
below 0.1—. Due to their similar form, Hamiltonian matrix and overlap matrix have the same 
eigenvectors \u^) = (^/2:^^ ^^2)' '^i^'^ V defined as f (k) = \f {k)\ e'"^^^^ = w {k) e'''^^^\ 
This yields the electronic eigenvalues 

ik) = ^i^. (19) 

We note that the overlap makes the energy bands asymmetric with respect to the Fermi 
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FIG. 6: (a) Optical matrix elements of the tight-binding model with overlap for the two- 
dimensional graphene. Since we are interested in the relative change of the momentum matrix 
elements with respect to overlap, arbitrary units are used, (b) Energy bands of the tight-binding 
model with overlap for the two-dimensional graphene. The values used are the following: 70 = 2.7 
eV, Ep = -5 eV, and sq =0.1. 

level and has large influence on bands. The full form of the Hamiltonian with overlap, H, is 

Ep - so'yow'^ (k) {'Jo - SoEp) f (k) 



H{k) 



(20) 



1 - slw^ (k) y (^^ _ ,^Ep) r (k) Ep - soiow' (k) 
Since Sq is less than 0.1 we can safely discard the prefactor in Eq. ()20|). The diagonal part 
of the Hamiltonian matrix is proportional to the unit matrix in both cases, with or without 
overlap, and it does not contribute to the interband momentum matrix element. Therefore, 
one can easily calculate the intraband momentum matrix element in a compact form as 



P 



- (m+I V// {k)\u) = — (70 - ^pso) 



w (k) 



-Vifi (k) 



(21) 



1 — Sgiw^ (A;) 

Eq. (|21|) shows us that for the most used parametrization {Ep = eV) , the overlap does 
not play any role on the inter-band momentum since sq is less than 0.1 and we can safely 
discard the second order term is Sq. However, in order to fit the experimental dielectric 
function with the nearest-neighbor model, orbital overlapping is invoked in Ref. It was 
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found that 70 = 2.7 eV and Ep = -5 eV by assuming sq =0.1. With this parameterization 
the numerical results are shown in Fig. IHl for the electronic bands and momentum matrix 
elements. The momentum matrix elements for the case with overlap are practically the 
same as those of the first order approximation Hamiltonian (Eq. (jll|) ). In the same time, 
the electronic bands generated by the first order Hamiltonian are different from those of the 
full Hamiltonian with overlap. 

IV. CONCLUSIONS 

We investigated the influence of the non-orthogonal orbitals on optical matrix elements in 
tight-binding models. A diagonal coordinate operator in the orthogonalized basis not only 
ensures the gauge invariance but also induces intra-atomic contributions to the coordinate 
operator in the original (atomlike and non-orthogonal) basis. Moreover, the Hamiltonian 
matrix in the orthogonal basis is longer ranged than the Hamiltonian matrix in the initial 
non-orthogonal basis. As a consequence, one can justify the nearest-neighbor interaction 
of the spin-orbit couplingi^. It enables to describe the Dresselhaus term, which is not 
considered in the usual treatment of the spin-orbit coupling2S. 

Simple models are analyzed. The flrst model studied was the monoatomic linear chain 
with two orbitals per site as an approximation to the Kronig-Penney model. The model 
was also used in Ref. Q to show the role played by the intra-atomic matrix elements of 
the momentum operator. We found that, although the tight-binding model with overlap 
exhibits almost the same energy bands as the one with orthogonal orbitals, the optical 
matrix elements are closer to the exact matrix elements of Kronig-Penney model. The 
second model studied was the biatomic linear chain with one orbital per site. This case 
showed that optical matrix elements decrease with overlap increasing. We also analyzed 
the optical matrix elements of the tight-binding model for two-dimensional graphite at low 
energies (between n and n* electronic states). Optical matrix elements remain unchanged 
with respect to the overlap when the usual parametrization Ep= eV is adopted, while the 
bands change drastically. However, non-vanishing orbital overlapping and Ep= -5.0 eV are 
needed for better agreement with experimental data^^. 

In complete analogy with the above arguments, one can use a piecewise constant coor- 
dinate operator in the orthogonal representation of tight-binding linear muffin-tin orbitals 
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methoda^i to calculate optical spectra from ab-initio. The procedure will be faster because it 
will enable to calculate optical spectra directly from energy band calculations by employing 
the fast Fourier transformation and without directly evaluating the momentum operator. 
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